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ABSTRACT 

We investigate the 21-cm signature that may arise from the intergalactic medium 
(IGM) prior to the epoch of full reionization (z > 5). In scenarios in which the IGM 
is reionized by discrete sources of photoionizing radiation, the neutral gas which has 
not yet been engulfed by an H n region may easily be preheated to temperatures well 
above that of the cosmic background radiation (CBR), rendering the IGM invisible 
in absorption against the CBR. We identify three possible preheating mechanisms: 
(1) photoelectric heating by soft X-rays from QSOs, (2) photoelectric heating by soft 
X-rays from early galactic halos, and (3) resonant scattering of the continuum UV 
radiation from an early generation of stars. We find that bright quasars with only 
a small fraction of the observed comoving density at z ~ 4 will suffice to preheat 
the entire universe at z ^ 6. We also show that, in a Cold Dark Matter dominated 
cosmology, the thermal bremsstrahlung radiation associated with collapsing galactic 
mass halos (10 10 - 10 11 M Q ), may warm the IGM to ~ 100 K by z ~ 7. Alternatively, 
the equivalent of ~ 10% of the star formation rate density in the local universe, 
whether in isolated pregalactic stars, dwarf, or normal galaxies, would be capable of 
heating the entire IGM to a temperature above that of the CBR by Lya scattering in 
a small fraction of the Hubble time at z ~ 6. 

In the presence of a sufficiently strong ambient flux of Lya photons, the hyperfine 
transition in the warmed H I will be excited. A beam differencing experiment would 
detect a patchwork of emission, both in frequency and in angle across the sky. This 
patchwork could serve as a valuable tool for understanding the epoch, nature, and 
sources of the reionization of the universe, and their implications for cosmology. We 
demonstrate that isolated QSOs will produce detectable signals at meter wavelengths 
within their "spheres of influence" over which they warm the IGM. As a result of 
the redshifted 21-cm radiation emitted by warm H I bubbles, the spectrum of the 
radio extragalactic background will display frequency structure with velocity widths 
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up to 10, 000 km s . Broad beam observations would reveal corresponding angular 
fluctuations in the sky intensity with 5T/T < 10~ 3 on scales 6 ~ 1°. This scale is 
set either by the "thermalization distance" from a QSO within which Lya pumping 
determines the spin temperature of the IGM or by the quasar lifetime. 

Radio measurements near 235 and 150 MHz, as will be possible in the near future 
using the Giant Metrewave Radio Telescope, may provide the first detection of a 
neutral IGM at 5 < z < 10. A next generation facility like the Square Kilometer Array 
Interferometer could effectively open much of the universe to a direct study of the 
reheating epoch, and possibly probe the transition from a neutral universe to one that 
is fully ionized. 

Subject headings: cosmology: theory - diffuse radiation - intergalactic medium - 
quasars: general 
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1. Introduction 

At epochs corresponding to z ~ 1000, the intergalactic medium (IGM) is expected to 
recombine and remain neutral until sources of radiation develop that are capable of reionizing 
it. The application of the Gunn-Peterson (1965) constraint on the amount of cosmologically 
distributed neutral hydrogen to QSO absorption spectra requires the universe to have been highly 
ionized by z ~ 5 (Schneider, Schmidt, & Gunn 1991). It thus appears that substantial sources 
of ionization were present at these or earlier epochs, perhaps QSOs and young galaxies, or some 
as yet undiscovered class of objects such as Pop III stars or decaying particles. Establishing 
the epoch of reionization is crucial for determining the impact of reionization on several key 
cosmological issues, from the role reionization plays in allowing collapsed objects to cool and 
make stars, to determining the small-scale structure in the temperature fluctuations of the cosmic 
background radiation (CBR) (e.g., Sugiyama, Silk, Sz Vittorio 1993; Bond 1995; Hu & White 
1996). Conversely, probing the epoch of reionization may provide a means of detecting the onset 
of the first generation of stars. 

Very little is known about the nature of the first bound objects and of the thermal state 
of the universe at early epochs. In particular, the history of the transition from a neutral IGM 
to an almost fully ionized one can reveal the character of the ionizing sources. If these are 
uniformly distributed, like an abundant population of pregalactic stars or decaying dark matter, 
the ionization and thermal state of the IGM will be the same everywhere at any given epoch, 
with the neutral fraction decreasing gradually with cosmic time. In a CDM dominated cosmology, 
this may be the case since bound objects sufficiently massive (~ 1O 6 M ) to make stars form at 
high redshift. If a sufficient fraction of their ionizing photons are able to escape, the stars in 
these objects would reionize the universe by z > 20 (Couchman & Rees 1986). In inhomogeneous 
reionization scenarios instead, widely separated but very luminous sources of photoionizing 
radiation such as QSOs, present at the time the IGM is largely neutral, will generate expanding 
H II regions whose size is only limited by the individual source lifetime. The universe will be 
divided into an ionized phase whose filling factor q{z) increases with time, and an ever shrinking 
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neutral phase (Arons & Wingert 1972; Meiksin & Madau 1993). If the ionizing sources are 
randomly distributed, the H n regions will be spatially isolated for q <C 1. The IGM will become 
completely reionized at the breakthrough epoch, when q = 1 and the IGM is fully transparent at 
the Lyman continuum. Current QSO counts suggest that the breakthrough epoch may occur as 
recently as z > 5 (Meiksin & Madau 1993). 

It is therefore of great interest, in testing rival reionization scenarios, to seek any evidence 
of non-uniformity in the ionization structure of the IGM by studying the thermal and ionization 
history of the intergalactic gas surrounding isolated sources of radiation (see, e.g., Miralda-Escude 
& Rees 1994). In this paper we shall discuss how radio-astronomical measurements of spectral 
and spatial structures in the 21-cm line emission and absorption from H I at high redshifts may 
provide a valuable tool for probing the epoch of the likely heat input at 5 ^ z ^ 10 and its nature, 
e.g., whether it was limited to small isolated regions or was relatively diffuse. We will confine our 
analysis to epochs prior to complete reionization, and address the physical state of the gas which 
has not yet been engulfed by an ionized region. 

For intergalactic H I to be observable in 21-cm line emission or absorption against the 
cosmic background radiation, the spin temperature must differ from the CBR temperature. As 
the former will in general be a weighted mean between the matter and CBR temperatures, an 
efficient mechanism for coupling the populations of the hyperfine levels to the kinetic state of the 
intergalactic gas is needed to unlock the spin state from the CBR. At the low particle densities 
characteristic of the IGM, the dominant coupling mechanism is level-mixing via Lya photon 
scattering - the "Wouthuysen-Field" effect. We present in this paper a new heating mechanism, 
arising from atomic recoil in response to the scattering of resonant photons. We demonstrate that 
because of resonant scattering heating, Lya coupling makes the detection of the IGM in absorption 
against the CBR quite unlikely, since the same Lya photons that pump the hyperfine levels heat 
the atoms as well. This is not entirely surprising, as the energetic demand for raising the kinetic 
temperature above that of the CBR is relatively small, ~ lCT 7 eVcm~ 3 at z ~ 7, hence even a 
relatively inefficient heating mechanism like Lya resonant scattering may become important. 
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We find that the heating requirements are quite generally met in a variety of scenarios 
in which radiation sources turn on at high redshifts. We show that, in models based on QSO 
photoionization, and well before the H n region network has fully percolated, the (mostly) neutral 
IGM between the H n bubbles will be photoelectrically heated to temperatures above a few 
hundred degrees by soft X-ray photons from the QSOs. In Cold Dark Matter (CDM) models a 
similar mechanism may efficiently warm the IGM at z < 7, since galactic mass halos with virial 
temperatures ~ 10 6 K are significant sources of thermal bremsstrahulung radiation. Alternatively, 
if protogalaxies at these early epochs are forming stars at a rate which is only ~ 10% of the 
global present-day star formation rate, heating by background Lya "continuum" photons could 
drive the kinetic temperature of the still neutral material well above the radiation temperature of 
the microwave background. As a consequence, we conclude that most of the neutral IGM will be 
available for detection at 21-cm only in emission. 

Two scenarios serve to illustrate the observability at meter wavelengths of the intergalactic 
medium during reionization. An isolated QSO turning on in a neutral medium will create an H II 
bubble surrounded by an X-ray warmed neutral zone. Sufficiently near the QSO, the spin state 
of the neutral IGM will be coupled to its kinetic state by redshifted continuum Lya photons from 
the QSO. This zone will then emit 21-cm radiation that would be detectable against the CBR as 
a differential intensity, in angle or frequency. A stronger signal may be achieved in the presence of 
a Lya continuum background sufficiently strong to couple the spin temperature of the cold and 
neutral IGM outside the QSO "sphere of influence" to the kinetic temperature of the IGM. The 
larger signal results from the enhanced contrast between the 21-cm absorbing medium far from 
the QSO, and the 21-cm emitting gas in the X-ray warmed neutral zone surrounding the QSO 
H II region. This phase, however, will be short-lived: through atomic recoil the Lya background 
photons will soon heat the IGM far from the QSO to a temperature above that of the CBR, and 
absorption will cease. 

The plan of the paper is as follows. In § 2 we discuss sources of radiation at high redshift. 
In § 3 we review the mechanisms that determine the level populations of the hyperfine state of 
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neutral hydrogen. In § 4 we treat the preheating of the IGM prior to complete reionization, and 
show that a 21-cm absorbing IGM can exist only as a transitory state. In § 5 we argue that the 
21-cm line signature from a neutral IGM, redshifted to meter wavelengths, will appear as a small 
feature in the radio extragalactic background on scales of < 1°, and may be detectable with the 
Giant Metrewave Radio Telescope. Finally, in § 6, we summarize our conclusions. 



In the following, we assume that the universe is completely reionized by discrete sources at 
z > 5, and that the IGM at such early epochs is uniform on large scales. Unless otherwise stated, 
we shall adopt a flat cosmology with qo = 0.5 and Hq = 50 /150 kms" 1 Mpc -1 . 



When an isolated point source of ionizing radiation turns on, the ionized volume initially 
grows in size at a rate fixed by the emission of UV photons, and an ionization front separating 
the ionized and neutral regions propagates into the neutral gas. The evolution of an expanding 
cosmological I-front in a uniform IGM is governed by the equation 



(Shapiro 1986), where rj is the proper radius of the I-front, «h is the hydrogen density of the IGM, 
H is the Hubble constant, S is the number of ionizing photons emitted by the central source per 
unit time, as — 2.6 x 10 _13 T 4 7°' 845 cm 3 s _1 is the recombination coefficient to the excited states of 
hydrogen, and T4 is the gas temperature in units of 10 4 K. A helium to hydrogen cosmic abundance 
ratio x = 1/12 is adopted, with both H and He fully ionized inside the H 11 region. The right 
hand side of equation (|l|) is a measure of the net flux of ionizing photons reaching the I-front after 
subtracting recombinations, and vanishes at the Stromgren radius, rst = {3S , /[47rn^as(l-|-x)]} 1 / 3 . 
Across the I-front the degree of ionization changes sharply on a distance of the order of the mean 



2. Sources of Radiation at High Redshift 



2.1. QSO H 11 Regions 




(1) 
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free path of an ionizing photon. Denoting by ^igm the baryonic mass density parameter of the 
intergalactic medium, the I-front will expand at the speed of light 

I S \ 1/2 , ,/ S \'' 2 /SliauhLY'' 2 /l + 2 \- 3 '2 

ri< (w) " (2 - 85Mpc) (i5*f0 pS**) (— ) ' (2) 

to subsequently slow down because of geometrical dilution. In an expanding IGM, the ionization 
front generated by a QSO, with S ~ 10 s _1 , or a star-forming galaxy, with S ~ 10 53 s _1 , will 
fail to grow to even half of its Stromgren radius (Shapiro 1986; Madau & Meiksin 1991). In this 
case radiative recombinations within the H n region may be neglected on the right-hand side 
of equation (|l|) (proportional to r| t — rf), to obtain at a short time after the turn-on epoch z 
(corresponding to the redshift interval Az -C 1 + z), 

and 



10 57 s -iy &u ^ Q 05 y ^ 7 
For comparison, the radius of the light front around each source is 



1 + z \ -5/2 

cAt ps (46.3 Mpc)/^ 1 ( —I- ) Az. (5) 



2.2. Virialized Halos 

In CDM dominated cosmologies, galactic mass halos will begin to collapse at high redshift 
(z > 5). The hot gas in these halos will produce a soft X-ray background radiation field via 
thermal bremsstrahlung. We estimate the number of halos using the Press-Schechter formalism 
(Press & Schechter 1974), which should suffice for an order of magnitude estimate for systems on 
these scales. 

The (proper) volume emissivity per unit frequency due to thermal bremsstrahlung radiation 
from isothermal, optically thin halos is 

e h ak» = (3.2 x lO^ergsMpc^s^Hz- 1 ) / dr N(r , z)(l + ^) 3 < (r )f r^ro)- 1 / 2 

Jo 

xexp[-hu/kT(r )], (6) 
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where N(tq,z) is the number of halos per unit comoving volume per initial comoving size rp at 
redshift z, nh a io( r o) is the internal hydrogen density of a halo, and T(r ) is its temperature. The 
hydrogen is assumed to be shock heated to the virial temperature and collisionally ionized. In 
spherical top-hat collapse, the temperature T and radius r v of a virialized halo are related to 
ro by kT(ro) = 1.4/jmH(l + z^HqTq and ro = 5.6(1 + z c )r v , where z c is the collapse redshift 
of the perturbation, and [i ~ 0.6 is the mean molecular weight of a fully ionized gas of cosmic 
abundances. The internal density of the cloud exceeds the external density at the projected 
time of collapse by a factor of 178. Normalizing to a baryon density of ^igm^Io = 0.05 gives 
rihaio = 1-8 x 1CP 5 (1 + z c ) 3 cm~ 3 . According to the Press-Schechter formalism, 

1.686 2 (l + z) 2 



1.686(1 + z) rflogA(r ) 
N{r ,z)dr = -3 — — — -—- exp 
(27r) 3 / 2 r^A(r ) dlog(r ) 



dr . (7) 



2A(r ) 2 

Here, A(ro) is the rms mass fluctuation in a sphere of radius ro. We will adopt a pure CDM 
spectrum for aqo = 0.5, Hq = 50kms _1 Mpc -1 universe (Bardeen et al. 1986). This model appears 
to provide an adequate representation of the power spectrum on small scales when appropriately 
renormalized, even though a higher normalization is needed for an optimal match to the 
temperature fluctuations measured by COBE on large scales (Gorski et al. 1995). Normalizing to 
the abundance of clusters, A(16Mpc) ps 0.6 (White et al. 1993), we find that for 0.2 < r < 1 Mpc, 
A(ro) ~ 3.6r , where ro is measured in megaparsecs. 

The total emissivity above a minimum photon energy hv m \ n is ehaio = [kT(ro)/h} 

^halo (^min) i 

after integrating equation (^) over frequency. While the emissivity from an individual halo 
increases for increasing temperature (for kT < hu m i n ), the number of halos decreases rapidly. At 
any given redshift, the peak contribution to the integral in equation @ arises therefore from 
some intermediate size halo. We may thus approximate this integral by performing a second order 
Taylor expansion of the exponent about the halo size. Most of the bremsstrahlung radiation is 
found to be emitted by halos with an initial comoving radius 

4/H fl + z \ -12/H 



^0 



/ hv ■ \ 4/ii f \ + 2f\ ~ 11 11 
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The total emissivity is given by 



-3 „-l 



e ha i « (2.6 x 10 4d ergsMpc _d s 



S 3/8 (i +z 



o 



7 



7 r 



fcT(fp) 
20ryd 



V20rydJ 



-1/2 



x exp 



-r,.:$[i^ 



,3/4 1 n fl+A 1 / fez/ min \ f -2 
— J l,207ydj r 



(9) 



We shall use this expression in § 4.2.3 below to estimate the virialized halo contribution to the 
heating rate of the IGM. 



2.3. Early Stellar Populations 

An early generation of stars may provide a background radiation field of Lya photons that 
arises from the collective redshifted continua of the stars. The background specific intensity 
associated with a constant comoving emissivity from the present epoch to redshift z max , as seen 
by an observer at redshift z, is given by 

J a (z) = - £ K,0)(l + z) 3 ^ / (i + z ')-<*s dz > (10) 

47r Jz dz 

where e(u a , 0) is the Lya volume emissivity due to stars in galaxies at the present epoch, dijdz' 
is the proper differential line element in a Friedmann cosmology, and we have assumed a power 
law energy spectrum. Here, 1 + z max ~ 4(1 + z)/3, as photons emitted at higher redshifts will 
have wavelength below 912 A at the source, and will be strongly attenuated by Lyman-continuum 
intergalactic absorption. Weakly dependent on as, equation (|l^) may be rewritten as 

e(v a ,0) 



J a [z) ss (nr^ergscm^s^Hz^sr- 1 ) 



10 25 ergs Mpc 3 s _1 Hz -1 



1 _|_ z \ 3/2 

-¥) ■ <»> 



Using the results of the Ha survey of Gallego et al. (1995), we estimate the continuum 
emissivity from present day galaxies at the frequency of Lya to be, to within a factor of 2, 
e(v a ,0) ~ 4 x 10 25 ergs Mpc~ 3 s _1 Hz -1 . Only a fraction of this emissivity would be available at 
high redshift from dwarf galaxies or the progenitors of today's bright galaxies. The Lya continuum 
emissivity at early epochs, however, may be even larger if pregalactic Pop III stars are present at 
these times. 
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The IGM will still be mostly neutral if the background intensity at the Lyman edge, Jl, 
satisfies 

Jl < ^n H (0)(l + z) 3 » (lO-^ergsem^s^Hz^sr- 1 ) {^f^j (^)\ ( 12 ) 

i.e., if less than one ionizing photon per hydrogen atom has reached the IGM by that epoch. 
Because the integrated UV spectrum of a stellar population is characterized by a strong 
Lyman-continuum break (about a factor of 4), and, more importantly, only a small fraction of 
ionizing photons might escape from collapsed star-forming regions into intergalactic space (cf. 
Madau &; Shull 1996), one expects Jl <C J a - Hence, based on a comparison of equations (p"l|) 
and (^), it is plausible to assume a largely neutral medium even in the presence of a significant 
Lya background from an early generation of stars, J a ,-2i ~ 1- 

3. Spin Temperature of the IGM 

The amount of emission or absorption from a neutral intergalactic medium is determined by 
the spin state of the atoms. The spin temperature may be coupled to the matter both through 
the scattering of Lya photons - the Wouthuysen-Field effect (Wouthuysen 1952; Field 1958) - 
and atomic collisions. At the low particle densities characteristic of intergalactic gas, and in the 
presence of radio-quiet radiation sources, the dominant coupling mechanism is Lya scattering. 

3.1. Lya Radiation Color Temperature 

The Wouthuysen-Field effect mixes the hyperfine levels of neutral hydrogen via the 
intermediate step of transitions to the 2p state. According to this mechanism, an atom initially 
in the n = 1 05*1/2 singlet state may absorb a Lya photon that will put it in the n = 2 \P\ji or 
1P3/2 state, allowing it to return to the triplet n = 1 1S1/2 state by spontaneous decay. Similarly, 
a triplet atom may absorb a slightly lower frequency photon to reach the same upper state, 
followed by decay to the singlet state. The efficiency of the effect is determined by the frequency 
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dependence of the background radiation field J u near Lya . In particular, the excitation rates 
of the hyperfine levels depend on the relative intensities of the blue and the red wings of Lya . 
To observe the 21-cm line in absorption or emission, the excitation and de-excitation rates of 
the hyperfine levels must be comparable to or exceed those due to 21-cm continuum background 
radiation. 

The Lya line absorption cross section is 

ire 2 3 

a v = a a <\>{y) = f a< j>( v ) = —\ 2 a A a cP(v) (13) 

(neglecting the very small correction due to stimulated emission, of order X^J a /2hc <C 1), 
where A a and f a are the spontaneous Einstein coefficient and upward oscillator strength for the 
transition, (f){v) is the line profile function (with normalization / (j)(y)dv = 1), and all other symbols 
have their usual meanings. The absorption rate of Lya photons is proportional to Bjk J du(f)(u)J v , 
where Bjk is the Einstein stimulated absorption coefficient from level j to level k, where we now 
mean to include the hyperfine structure of the atom. [The coefficients are related to the total 
spontaneous decay rates A^j between the various hyperfine levels by Bjk = (\^ k /8irh)(gk/ 'gj)A^\ 
When the intensity varies smoothly across the H I absorption profile for Lya , that is, over at least 
one Doppler width Az/q = (b/c)v a , where b = (2kTx /ma) 1 / 2 and v a is the line center frequency, 
it is useful to express the effect of the mixing mechanism in terms of the color temperature T a of 
the radiation field, which we define by 

kT a dhv ' 1 ' 

evaluated at v = v a , where M(y) is the photon occupation number at frequency v, 
N{v) = c 2 J u / (2hv 3 ) . If the background continuum spectrum near Lya is a power law, 
J u oc v~ a s, then T a = (hv a /k)/(3 + a s ) ~ 3 x 10 4 K[4/(3 + a s )]. It has been shown by Field 
(1959b), however, that, because of the large cross section for resonant scattering, the shape of 
the radiation spectrum near Lya will be determined in all the relevant physical situations by the 
kinetic temperature of the scattering material. This tendency for thermalization is caused by the 
recoil of the atom, which shifts the profile to the red. As stimulated emission is negligible, the 
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radiation behaves like a classical relativistic Boltzmann gas, and the specific number density of 
photons, n v = A7rJ u /chu, assumes the quasi-LTE form n v = const x exp[— h(y — iy a )/kTx], where 
Tk is the kinetic temperature of the gas, to which the color temperature relaxes. 

The influence of cosmological expansion on the Lya profile has been discussed by Field 
(1959a). The effect of expansion is determined by two timescales, the mean- free scattering time 
for a Lya photon, b/ (cnm&a^a)-, where nm is the density of neutral hydrogen, and the time it 
takes the Hubble expansion to redshift a photon through the resonance, b/(Hc). Their ratio, also 
known within the framework of line radiation transfer theory as the Sobolev parameter (and equal 
in this context to the inverse of the Lya optical depth), is 



As this is much smaller than unity, the Lya part of the radiation background spectrum establishes 
close thermal contact with the expanding IGM, and a quasi-static treatment of the cosmological 
line transfer is an adequate approximation (see, e.g., Rybicki & Dell' Antonio 1994). The profile of 
Lya radiation in an expanding cloud has been calculated in detail by Deguchi & Watson (1985). 
These authors consider the case in which Lya photons are created by recombinations inside the 
cloud. They show that T a — > Tk for 7 < 10~ 5 , similar to what is found in the absence of velocity 
gradients (Field 1959b). The color temperature approaches the gas kinetic temperature even more 
quickly if the cloud is irradiated by continuum radiation in the vicinity of Lya , as this spectrum 
is already flat in frequency, and only a small number of scatterings is needed for the radiation field 
to achieve the proper slope (Deguchi & Watson 1985). In the next section we shall see how, in the 
presence of a strong ambient flux of Lya photons, the color temperature governs the hyperfine 
populations of atomic hydrogen. 

3.2. Lya Photon Pumping 

In the absence of collisions, the spin temperature, T5, of the hyperfine levels of intergalactic 
hydrogen (as defined by the Boltzmann equation) is determined solely by the scattering of ambient 




(15) 
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Lya photons and by the microwave background radiation field at 21-cm (Field 1958, 1959a). Q In 
a steady state we have 

T s = r— , (16) 

1 + Va 

where Tcbr = 2.73(1 + z) K is the temperature of the cosmic background radiation (Mather et al. 
1994), and 

Voi = -J—7f- (17) 

^10 l-ot 

is the normalized probability, or Lya pumping "efficiency". Here, = hvi$/k = 0.07 K, 
A\q = 2.9 x 10~ 15 s" 1 is the spontaneous decay rate of the hyperfine transition of atomic hydrogen, 
Pio is the indirect de-excitation rate of the triplet via absorption of a Lya photon to the n = 2 
level, and we have assumed T$ ^> T*. In the absence of Lya pumping the spin temperature goes to 
equilibrium with the 21-cm background radiation field on a timescale T*/(Tcbr^4io) ~5x 10 4 yr, 
and neutral intergalactic hydrogen cannot be seen in either 21-cm emission or absorption. If y a 
is large, Tg — > T a ~ Tx, signifying equilibrium with the matter. The 1—^0 transition rate via 
Lya can be replaced by the total rate P a , 

P a = ( dn [ -?-(J v dv, (18) 
J J hv 

at which Lya photons are scattered by an H atom in the gas, Pio = 4P a /27. For a point source, 
I v = f u 5(£l), where f u is the specific flux from the source. For an isotropic radiation field, I v = J u . 
In the limit T a 3> TcbRj the fractional deviation in a steady state of the spin temperature from 
the temperature of the microwave background is 



Ts — Tqbr 



1 



TcBR 



-i -1 



DaTa 



(19) 



T s 

There exists then a critical value of P a which, if much exceeded, would drive T$ — > T a . This 
"thermalization rate" is found to be 

Ah - *™g?t** « (5.3 x 10-- (1±£) . (20) 



2 This is not strictly true in the presence of bright continuum sources of 21-cm radiation, such 
as radio-loud quasars (see § 3.5.). 
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When P a = P t h, we have T$ = 2Tcbr/(1 + TcBn/T a ). It should be emphasized that, while 
T a depends on the spectral gradient across line center, y a T a , and hence P t h> are only functions 
of the magnitude of J a , independent of the uncertainties regarding the actual profile of the 
Lya background photon spectrum. 



There are four possible sources of background Lya photons propagating into the (mostly) 
neutral material which has not yet been engulfed by an H n region: a) UV, non-ionizing 
continuum emitted by the radiating sources and then redshifted by the Hubble expansion to the 
Lya wavelength in the absorbing medium rest-frame; b) Lya photons generated by radiative 
recombinations and escaping from the H n regions; c) Lya produced locally in the warm neutral 
material surrounding the ionized gas by radiative recombinations; and d) by collisional excitations. 
Of these, only the first, continuum photons redshifted into local Lya , is a significant source of 
coupling, and we treat it next. We show that the remaining three are generally negligible in 
Appendix A. 

In the limit where the scattering occurs close to the emission redshift, we find for a point 
source 



where L at 4j is the source Lya luminosity per octave in photon energy in units of 10 ergs s , 
and rMpc is the (physical) distance in Mpc between emission and scattering. Thus, within several 
Mpc from a bright, radio-quiet QSO, P a > P t h and the spin temperature of neutral hydrogen will 
be determined by the Lya continuum photons emitted from the quasar itself. 

Alternatively, we may consider the effect of an integrated intensity from a distribution of 
sources. The Lya scattering rate is then 



3.3. Level Mixing by Lya Continuum Photons 



P a « (2.3 x 1(T 10 s- 1 ) r M 2 L aA7 , 



(21) 



P a « (8.5 x 1(T 12 s- 1 ) J a _2i, 



(22) 



where J a ,-2i is the intensity of the background diffuse field at Lya in units of 
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10~ 21 ergs cm -2 s _1 Hz -1 sr . For J Qj _2i 1, the spin system will go into equilibrium 
with the microwave background. Equation ( p2[ ) is also approximately valid in the case of 
recombination Lya photons, as radiative transfer calculations predict dJ v /dv = at the Doppler 
core of the transition (Urbaniak & Wolfe 1981). Spectral gradients induced by atomic recoil only 
result in higher order corrections for P a . The rate is therefore dependent only on J a . 

The Lya pumping mechanism poses a difficulty for detecting the IGM in absorption against 
the CBR, since, as we will show below, the same Lya photons that mix the hyperfine levels of 
hydrogen heat the atoms as well. Before we discuss this problem, we briefly address the influence 
of particle collisions and direct radio continuum on the spin temperature. 



3.4. Hyperfine Excitation via Particle Collisions 

We demonstrate in this section that, except at very high redshifts, or in overdense regions, 
collisions are not effective in coupling the spin temperature to the kinetic temperature. One can 
show that, in the general case, the spin temperature may be written in the form (Field 1958) 

rp s _ Tcbr + UaTg + VcTfC ^23) 
l+Va + Vc ' 



where 



(24) 

^10 J-K 



and Cio is the rate of collisional de-excitation of the triplet level (cf. eqs. 16] and |17]]). Emission 
occurs for 

UaT a + UcTk m - m /oC\ 

■ *iT K > T C BR, (25) 

Va + Vc 

and absorption otherwise. Unless the electron density is larger than a few percent of the hydrogen 
density, the collisional de-excitation rate is dominated by collisions with other hydrogen atoms 
through spin-exchange (Purcell & Field 1956). Allison & Dalgarno (1969) give a maximum 
de-excitation rate of Cio/tih = 3.3 x 10~ 10 cm 3 s _1 for Tx < 1000 K, with a value a factor 2.6 
lower at T% = 100 K. To couple the spin temperature to the kinetic temperature requires a 
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collision rate C\q > -Aio^cbr/^*) or a density «h > 3 x 10 4 (1 + z) cm 3 . This may occur only 
for 1 + z > 55(^igm^5o/ -0 5 ) _1/2 - At z ~ 6, collisions will mix the hyperfine levels if Sp/p > 100, 
i.e., in the case of nearly virialized halos. 



3.5. Radio Loud Quasars 

So far, we have limited our attention to radio quiet sources of radiation. However, 21-cm 
radiation emitted by quasi-stellar radio sources could raise the spin temperature of neighboring 
intergalactic H I above the CBR temperature (Bahcall & Ekers 1969; Urbaniak &: Wolfe 1981). At 
a distance r from the QSO, the brightness temperature at 21-cm is 

c 2 L 

TR= 8^f 4^> (26) 
where L u is the absolute luminosity of the radio-emitting QSO at 1.4 GHz. The typical radio 
quasar power is a small fraction of the B-band luminosity, 

(vL u ) 5GUz ^4xW- 4 (vL u ) B (27) 

(e.g., Phinney 1985). So, assuming vL v is roughly constant in the radio, we derive 

{yL v ) B 



T R «(30K) 



10 47 ergs s 



7-Mpc- (28) 



This temperature must be added to Tcbr in equation fllq ) for determining the level populations. 
Comparison with equation (^T|), however, shows that in the presence of UV photons this term is 
very small compared to the Lya pumping term and so can safely be neglected. 



4. Preheating of the IGM 

While a cold IGM may be detected in absorption against the CBR, when the spin and 
kinetic temperatures are coupled through the Wouthuysen-Field mechanism, the period for which 
absorption is possible is short-lived. The reason is that the Lya photons that couple the spin and 
kinetic temperatures will quickly heat the IGM through atomic recoil to a temperature above that 
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of the CBR. In this section, we return to the radiation sources introduced in § 2 and discuss their 
roles as heating sources by Lya scattering and by photoelectric heating by soft X-rays. 



4.1. Lya Heating 

The effect of recoil due to the finite mass of the atom was first discussed in resonance-line 
scattering by Field (1959b) (see also Adams 1971; Basko 1981). Contrary to the case of radiation 
scattered by free electrons (where a systematic shift in frequency due to recoil can produce large 
cumulative effects) , in problems of resonance- line transfer a photon on the red side of the line has 
a heavy bias toward scattering back to the blue. This bias suppresses the additive effect of recoil: 
because of Doppler redistribution, a photon largely loses memory of its injection energy, and the 
background intensity develops only a slight asymmetry in frequency at line center (matching the 
slope of the Planck function in the Wien limit, see Field 1959b) even at the large Lya optical 
depths normally encountered in astrophysical situations. However, it is through the recoil of the 
atom that the radiation field exchanges energy with the gas in scattering problems. In this section 
we shall discuss in some detail the effect of repeated Lya resonant scattering on the thermal 
state of a cold, neutral IGM. A more formal derivation of the rate of transfer of energy from the 
radiation field to the atoms is provided in Appendix B. 

The average relative change in a Lya photon's energy E after having been scattered by a 
hydrogen atom at rest is 

O = -^V-10" 8 , (29) 
E nine 1 

where ran is the mass of the hydrogen atom. It should be noted that this is an approximation 
valid only for hu a » kTx- In the opposite limit, energy will flow from the atoms to the photons. 
Through recoil, energy is transferred from photons to atoms at a rate 

AE 

E a = -{—)hv a P a . (30) 
where P a is the Lya scattering rate per H atom. In the case of excitation at the thermalization 
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rate P t h, equation (|30|) becomes 



^ _ K (220 K Gyr-) (l±i) . (31) 

This rate is sufficient to drive the intergalactic gas kinetic temperature above the temperature of the 
cosmic background radiation in a fraction ~ 15% of the Hubble time at z = 6. As a consequence, a 
low density IGM is observable in absorption at 21-cm for only a brief interval of time once coupling 
is established. In general, coupling through the Wouthuysen-Field mechanism leads to emission. 

This may be illustrated for the case of an early population of stars. It follows from 
equations (11) and (22) that the equivalent of only ~ 10% of the global present-day star formation 



rate is required for P a to exceed P t h, and so drive T$ — ► Tr-. The IGM would initially be observed 
in absorption, but after ~ 10 8 yr, the IGM will be heated to a temperature above Tcbr by resonant 
scattering of the Lya photons. 

It is of interest to compare Lya heating with heating by X-ray photons. The amount of 
energy deposited via photoionization by soft X-rays is, allowing for secondary ionizations, a few 
tens of eV per photon. There are also a few free electrons for every photon, as recombinations 
are negligible. In contrast, the total energy available for Lya heating (whether the Lya photon 
is part of the continuum or due to recombinations) is, neglecting the energy lost to cosmological 
expansion, given by the photon energy multiplied by the fractional line width. This is typically 
0.001 eV per photon. Therefore, at least of order one Lya photon per atom is needed to raise the 
kinetic temperature of the IGM above Tcbr- If we denote by Lx the soft X-ray luminosity of a 
young galaxy (associated with, e.g., supernova remnants, X-ray binaries, etc.), then photoelectric 
heating will dominate over Lya heating if Lx/L a > 10~ 4 . As the typical observed value for L* 
galaxies is Lx/L a ~ 10~ 2 (e.g., Fabbiano 1989; Bennet et al. 1994), Lya photons from galaxies 
are an inefficient heating source compared to X-rays. ^| 



3 As shown by Rees (1985), the radiation pressure associated with trapped Lya background 
photons is also unlikely to affect the condensing out of overdense regions of the IGM even in the 
linear regime. 
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4.2. Soft X-ray Heating 

4-2.1. Heating Rate by Quasars 

It is widely believed that the integrated ultraviolet flux arising from QSOs and/or hot, 
massive stars in metal-producing young galaxies is responsible for maintaining the intergalactic 
diffuse gas and the Lya forest clouds in a highly ionized state at z < 5. In this and the next section 
we shall focus on a reionization scenario dominated by hard-spectrum quasar sources. In this case, 
while the radiation just shortward of the hydrogen Lyman edge will be absorbed at the ionization 
front generated by an individual QSO, very short wavelength photons will be able to propagate to 
much greater distances, r c > rj. The intensity of the QSO radiation flux at the light front will be 
reduced by the opacity of the (mostly) neutral gas surrounding the expanding H 11 region, 

Tv = n m (r c -r I )(o*? + xa **y ( 32 ) 

where cr™ and a^ el are the hydrogen and helium photoionization cross sections. For 
^igm^Io = 0-05, z > 5, and r c — rj = lOMpc, the photoelectric optical depth will drop below 
unity for photon energies above 16 ryd. Thus, most of the photoelectric heating of the IGM will 
be accomplished by soft X-rays. The heating by such high energy photons, however, is not as 
efficient as for photons near the photoelectric threshold because of losses to collisional excitation of 
H 1 and He 1 and to the collisional production of secondary electrons (Spitzer & Scott 1969). The 
fraction / of photon energy converted into heat is then a sensitive function of the H 11 fraction. At 
the light front, this fraction is the residual ionization left over from the recombination epoch. For 
h 50 = 1 and Oigm = 0.05, this is ~ 5 x 10~ 4 (Peebles 1993), which gives / « 0.17 (Shull k van 
Steenberg 1985). For a quasar with a spectral index as shortward of lryd, the heating rate per 
hydrogen atom at the light front is given by 

S r°° 

Ex(r c ) = a s^hu L f ^ dxx-'-^ia^ix - 1) + X ct? cI (x - \.%)\e~ T % (33) 

where x = vjv^, and vl is the frequency at the H 1 Lyman edge. Both the H 1 and 
He 1 photoelectric cross sections scale approximately as v~ 3 far from threshold. In this 
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approximation the integral becomes exact, and may be expressed as 

*<W « [(1400 KGn-'wr (^) rj-""] (A) („frr) (W^'" 

x(^)- 2 - <,S (l--)- (2+<,S)/3 ( 5 rffe)- ,8+ ° S,/S , (34) 

where T(x) is the usual gamma function. Here, r* normalizes the opacity of the IGM for soft 
X-ray frequencies according to r v = T*(fiiGM^5o/0-05)[(l + ,z)/7] 3 (r c /25Mpc)(l — rj/r c )(u/uL)~ 3 . 
Numerically, r* ~ 7.8 x 10 4 . For as = 1.5, the coefficient in square brackets is ~ 300 KGyr -1 . 
The heating rate is dominated by helium absorption, which exceeds the hydrogen contribution by 
~ 2:1. The intergalactic gas at the light front will be heated to a temperature above that of the 
cosmic background radiation in a fraction of the Hubble time, 

Athcat 



(a s -l)/3 
0.02^ 



(2+a S )/3 

_ /,_.. I ! II " ' i" I 

a s r{(2+ as )/3} 

9/2+as / n(2+c* s )/3/ r \ (8+a s )/3 



ft 50^ .l 7 J (,10^ i=l J V 0.05 ) 



For as = 1.5, the coefficient in square brackets is ~ 0.1. The H n region produced by a QSO will 
therefore be preceded by a warming front. Heating by soft X-rays is so efficient that the size of the 
warm zone will extend nearly to the light front. Note that, as the X-ray-heated bubbles around 
QSOs will survive as fossils even after the quasar has died, several generations N g of quasars may 
actually be responsible for preheating the entire IGM. For a typical QSO age of tQ 3 x 10 8 iV 9 ^ 3 
yr, the required QSO comoving space density to heat the entire IGM to a temperature above that 
of the CBR by z ~ 6 is ~ 10~ 10 Mpc^iV^ 1 . By comparison, the comoving space density of bright 
QSOs at z = 4 is ~ 100iV 9 times larger (Warren, Hewett, & Osmer 1994). Note that, if all bright 
galaxies undergo a quasar phase, QSOs must have a very short lifetime, and N g ~ 100. Soft X-rays 
from a few bright QSO sources could then prevent collapsing structures, such as protoclusters 
while still in the linear regime, from being detected in 21-cm absorption against the CBR (Scott 
& Rees 1990). 
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4-2.2. Numerical Solution for Expanding Warm Hi Regions Around Quasars 

We illustrate the thermal effect of an expanding H n region on the surrounding medium by 
numerically solving the time-dependent equations for photoionization and heating around a point 
source which turns on suddenly in a homogeneous, expanding universe. We have integrated as 
a function of time the ionization equations for hydrogen and helium, including photoionization, 
collisional ionization, and radiative recombination, on a radial grid with a local ionization flux 
equal to L v exp(— T u )/Ai:r 2 , where r v is given by equation (|32[). The flux is set to zero beyond the 
light front to prevent unphysical excess heating at large distances. We use the time-dependent 
ionization algorithm of Meiksin (1994), and include here the generation of secondary electrons 
by soft X-rays. We adopt the fitting formulas of Shull & van Steenberg (1985) to include the 
additional ionization and energy losses associated with the production of these electrons. We solve 
the time-dependent energy equation along with the ionization equations for a gas of primordial 
hydrogen and helium composition, including photoelectric heating and radiative losses due to 
collisional excitation and ionization, radiative recombination, Compton cooling off the CBR, and 
cooling by expansion. An accurate numerical solution for the growth of the H n region requires 
grid zones that are optically thin to the ionizing radiation when the hydrogen is still fully neutral. 
Thicker zones are ionized too rapidly, resulting in a too rapidly growing I- front. Because we wish 
to follow the heating out to the light radius, however, this restriction imposes a prohibitive demand 
on our computing resources. The demand imposed by helium ionization is less restrictive because 
of its lower density, and so may be met. We therefore adopt a compromise. A grid sufficiently 
fine to resolve the helium ionization is used, but it is somewhat too coarse to follow the evolution 
of the H II I- front accurately. To prevent too rapid a growth of the H n front, we set the flux to 
zero for energies less than 1.8 ryd, the He I ionization potential. We find this results in accurate 
placement of the H n front, though the details of its structure may not be quite correct. Because 
helium ionization and the soft X-ray heating are governed by the opacity outside the H n front, 
this approach yields an accurate solution for the He n and He in fronts and the external heating 
rate in the surrounding neutral IGM. 
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In Figure 1, we show the results of an integration for a QSO turning on at z = 6 with a 
gradual rise over a time of 10 7 yr to a peak photoionizing luminosity of S = 10 57 s _1 , and with 
as = 1.5, in a medium with fiiGM^io = 0.05. By z = 5.5, the QSO has heated the surrounding gas 
to a temperature exceeding that of the CBR to a (proper) distance of 23 Mpc, compared to the 
light front radius of 25 Mpc. The temperature climbs inward, reaching 1000 K at a distance of 10 
Mpc from the QSO. The ionization fronts are at 7.5 Mpc at this time. Neither the He n nor the 
He in fronts extend much beyond that of hydrogen, but for different reasons. Because of the much 
higher photon energy required to photoionize He n compared to H I , the radius of the He in front 
can exceed that of H n by at most a factor of (4~ as /x)^' 3 ~ 1.1, despite the low helium abundance 
(cf. Madau & Meiksin 1994). While the energy needed to photoionize He I is smaller, the IGM is 
more optically thick at the lower energy. At z = 5.5, the opacity of H I at the He I edge reaches 
unity a distance 10 kpc beyond the H n ionization front. Outside the H n region, H I efficiently 
shields the helium from the ionizing radiation above 1.8 ryd, ensuring that it remains mostly 
neutral. The He I is then able to warm through the absorption of the incident soft X-rays emitted 
by the QSO. 

4-2.3. Heating Rate by Virialized Halos 

Galactic mass halos provide another source of heating in CDM-dominated cosmologies. While 
before stars form the numbers of these halos are small, too small to account for the ionization 
of the IGM by UV photons, we demonstrate here that a sufficient number of halos will collapse 
to provide the soft X-rays required to heat the IGM to a temperature above that of the cosmic 
background radiation. 

We start with equation (0), which gives the background emissivity above a minimum photon 
energy. To specify its value, we observe that, at sufficiently low energies, the distance a photon 
travels before being absorbed in the IGM will be shorter than the mean separation of the sources, 
so that only part of the IGM will be penetrated by such photons. We set the minimum energy 
according to the criterion that the filling factor of photons with energy above /if m in be greater 
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than unity. For 5 < z < 10, this yields 



1 + z\ 3 / 2 

^ min «15.2ryd(^— , (36) 



7 

corresponding to a halo of comoving size ~ 0.9 Mpc, virial temperature ~ 1.5 x 10 6 K, and total 
mass ~ 2 x 10 11 Mq. Approximated in this manner, the emissivity is only half of the value given 
in equation (^), since halos with ro < ro are excluded by the filling factor criterion. Note that the 
collisional ionization level of helium in such hot halos is so large that absorption of soft X-rays 
within the halos themselves is completely negligible. Substituting equation ( |36| ) into equation (||), 
and taking into account the factor of 2, we find for the heating rate of the halos per hydrogen atom 



E)?° « (2.6 x lO^Gyr- 1 ) (X) ^±lj 'exp 



-6.6 



(37) 



s / 1 + z 

~1 

where / ~ 0.2 accounts for the energy loss to collisional excitations and ionizations per incident 
X-ray photon (Shull & van Steenberg 1985), for an assumed H n fraction of 10 -3 . We show the 
heating rate in Figure 2. The IGM will be heated to a temperature above that of the CBR in a 
fraction of the Hubble time, 



6.6 



1 + zV 



7 
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(7 .5xWxr7i±^ 

*H V0.27 V 7 

At z = 6, this fraction is 0.05. Because the merging of small halos into larger ones is an ongoing 
process, the full Hubble time will not be available to the halos at any given stage in the hierarchy 
for heating the IGM. Moreover, the halos will cool by radiative processes on a timescale of 
tcooi/tn ^ 0.7[(1 + z)/7]" 3/2 r 6 L8 (Gould & Thakur 1970), and even more quickly at high redshifts 
by Compton cooling off the CBR. The cooling, however, may not be significant for z < 10 before 
the gas is reheated by mergers in the next stage of the hierarchy. Allowing the halos to provide 
X-ray photons for 50% of the Hubble time, they will heat the IGM to a temperature of ~ 100 K 
by z = 6. At higher redshifts, halo heating may still be important in overdense environments like 
protoclusters where the number density of halos may be substantially boosted by statistical bias 
(Bardeen et al. 1986). Of course once stars start to form, they or, more likely, their remnants 
may dominate the soft X-ray emission of an individual halo. In this sense, the halo emissivity 
calculations outlined above present a lower limit to the heating rate of the IGM at early epochs. 
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It is worth noting that while the X-rays are an efficient heating mechanism of the gas, they 
are an inefficient source of photoionization. This is because the time required to photoionize the 
IGM by X-ray photons is larger by a factor of ~ (15 ryd/fcTcBR,) ~ 10 5 , since now it is the number 
of ionizing photons rather than their energy that counts. The ionization time thus exceeds the 
Hubble time by a few orders of magnitude, and the IGM will remain mostly neutral. 



4-2.4- Thermal Coupling Between Neutral Atoms and Ions 

In the discussion above, we have assumed that the energy deposited by the photons is strongly 
coupled to the hydrogen. The thermal energy, however, is initially contained in the photoejected 
electrons and ions only. In this section we shall verify that the equilibration mechanisms between 
the various ions and neutrals in the low density IGM are sufficiently rapid to ensure that all 
species are thermally coupled. 

We first consider Coulomb coupling between the ions. The time for two different species i and 
j to reach thermal equilibrium through Coulomb scattering is 

55 ' 5 - 87s) („, + 4^1ogA (l + ■ < 39 » 

(Spitzer 1962). Each species has a temperature T, density n, charge Z, and molecular weight A 
(A = 1/1836 for electrons). (This expression applies strictly only when the total energy between 
the two species, n{Fi + rijTj, is conserved during the equilibration process. For a single species, 
j = i may be set in equation (^) if the coefficient 5.87 is replaced by 8.06.) In the physical regime 
of interest, 20 < log A < 50. Once the electrons and protons each reach equilibrium separately 
(the shortest two timescales), the bottleneck for all the ions to relax to a single temperature is 
proton-electron scattering, with a coupling time 

Me.p) * (3.3 Myr) (^ (I±^ (^)V, (40) 

where x e is the number of electrons per hydrogen nucleus. After these two species equilibrate, 
singly and doubly ionized helium follow shortly through proton scatters. Note that, as the 
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post-photoionization (by 16 ryd soft X-rays) temperature of electrons is ~ 10 6 K, the time for 
protons to thermalize with such high energy electrons by Coulomb collisions alone actually exceeds 
the Hubble time. At these energies, however, the primary photoelectrons will quickly cool via the 
production of secondary electrons by the collisional ionization and excitation of neutral hydrogen 
(Spitzer & Scott 1969). The cooling rates of an electron of energy E ryd by these two processes 
are comparable, and are given by 



(Shull & van Steenberg 1985). For E = 16, this gives t coo \ ~ 0.6 Myr. The result is an energy 
distribution of primary and secondary electrons which declines sharply above the Lya excitation 
energy of 10.2 eV (Bergeron & Souffrin 1971). 

The coupling rate is still small for electrons of these energies, but they further lose energy 
through elastic scattering off neutral hydrogen atoms. The electron- hydrogen cross section for a 
10 eV electron is about 77toq, where ag is a Bohr radius (Moiseiwitsch 1962). The corresponding 
relaxation timescale is 



The electrons can finally cool to temperatures below ~ 10 4 K. Once a thermal distribution of 
electrons has been established, new energetic electrons generated by X-ray photoionization will 
reach equipartition with the cooler electrons on a much shorter timescale by Coulomb losses (Shull 
& van Steenberg 1985). 

The distance from the light-front for which the temperature of the neutral hydrogen atoms 
lags behind the electron temperature corresponding to equation (|42|) is ~ 8Mpc. This, however, 
assumes complete thermalization. The less stringent requirement that the neutral hydrogen atoms 
be heated to a temperature large compared to the cosmic background radiation temperature will 
be satisfied on a somewhat shorter timescale. The number of QSOs required to heat the IGM 
above the cosmic background radiation temperature may then exceed the estimate given in § 4.2.1 
by a factor of a few only. 




(41) 




(42) 
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The neutrals will reach a single common temperature relatively quickly. The interaction 
between neutral hydrogen atoms is dominated by elastic scattering. ^ The equilibration time is 
t eq = 2/nHx(ucr), where {ua) is the slowing down coefficient for hydrogen atom collisions (Spitzer 
1978). Using the elastic collision cross section from Allison & Smith (1971), this gives 

*eq(H, H) « (0.4 Myr) (±±f) " (^f^) (43) 

Thus all the species will reach a common temperature within a timescale much shorter than the 
Hubble time, and it is safe to assume that the hydrogen atoms may be described by a single 
Maxwellian energy distribution. 

5. Observations of Redshifted 21-cm Line Radiation 

In this section we discuss how radio measurements may reveal 21-cm line radiation from a 
neutral IGM at high redshift. The Giant Metrewave Radio Telescope (GMRT), scheduled for 
completion by 1997, has two frequency bands which are best suitable for detecting 21-cm radiation 
at z > 5, viz. 150 ± 20 and 235 ± 8 MHz (Swarup 1996). The central frequencies of these bands 
correspond to z = 8.5 and 5.0, respectively. Because of its high sensitivity for both compact and 
extended sources, the GMRT will be a valuable instrument for cosmological studies in the redshift 
ranges 4.85 < z < 5.25 and 7.35 < z < 9.9. We focus the discussion on measurements in the 150 
MHz band, corresponding to z ~ 8. 

We have assumed that the IGM at early epochs is uniform on large scales. One general 
issue is the effect of small-scale clumping. In a CDM-dominated cosmology, inhomogeneity 
on subgalactic scales develops before the first quasars form (if the latter require galactic-mass 
collapsed structures). This may delay the eventual percolation of H n regions, and increases the 
value of J needed to balance recombinations and maintain full ionization. It is unlikely that 

4 Provided that statistical equilibrium has been established between the singlet and triplet ground 
state levels. This is achieved on a timescale ~ 5 x 10 4 yr by collisions with CBR photons. 
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clumping would change the present discussion by very much. It leaves unchanged the amount 
of heating and ionization due to each UV photon, as well as the rate at which the heating and 
ionization fronts advance around a quasar. Even in an overdense region (unless the overdensity is 
extreme) , recombinations are unimportant over the lifetime of a single quasar (or during the time 
it takes for the heating front to advance) . It is also worth emphasizing that there will unavoidably 
be nonuniformities in the 21-cm intensity due to density inhomogeneities, but that these are 
restricted to smaller scales than those due to the patchy heat input (at least in the QSO case). 
While we confine our discussion to larger scales, it would nonetheless be worthwhile to consider 
means for their detection. 

To illustrate the basic principle of the observations we propose, consider a region of neutral 
material with spin temperature T$ 7^ Tcbr> having angular size on the sky which is large compared 
to a beamwidth, and radial velocity extent due to the Hubble expansion which is larger than the 
bandwidth. Its intergalactic optical depth at 21(1 + z) cm along the line of sight, 

3c 3 h 3 n m (0)A 10 15 ln - 2 .9,-l { TcBR \ AjlGM^jp A n q/ 2 (AA\ 

will typically be much less than unity. The experiment envisaged consists of two measurements, 
separated in either angle or frequency, such that one measurement, the fiducial, detects no line 
feature, either because there is no H 1 or because T$ « Tcbr, and the second at T$ / Tcbr- As 
the brightness temperature through the IGM is T& = TcBRe~ T + Ts(l — e~ r ), the differential 
antenna temperature (which is the relevant quantity since we are considering a detector operating 
as a comparison system), observed at Earth between this region and the CBR will be 



ST = (1 + z)-\T s - Tcbr)(1 - O « (0.011 K)h^ (^^) i^)^ ^' (45) 
where we have defined a "21-cm radiation efficiency" as 

^xhi(^^). (46) 

Here xhi refers to the neutral fraction of the hydrogen in the region for which T5 7^ ?cbr- As long 
as Tg is much larger than Tcbr (hence if there has been significant preheating of the intergalactic 
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gas), T) — ► xhi, and the IGM can be observed in emission at a level which is independent of 
the exact value of Tg. By contrast, when Tqbr 3> Is (negligible preheating), the differential 
antenna temperature appears, in absorption, a factor ~ Tqbr/Ts larger than in emission, and it 
becomes relatively easier to detect intergalactic neutral hydrogen (Scott & Rees 1990). Note that 
in a universe with go ^ 1/2, ST increases more nearly linearly with (1 + z), and the numerical 
coefficient in equation (|45|) may be larger by up to a factor 3[(1 + z)/9] 1//2 . Depending on the 
cosmology, 8T will increase toward low frequencies at a rate between v~ x l 2 and v~ x . 

The discrete nature of the sources of preheating and 21-cm excitation will cause 5T to have 
a distinctive structure both in frequency and in angle across the sky, as will be shown below. 
It is important to point out that the radio signal from the IGM will be swamped by the much 
stronger non-thermal background that dominates the radio sky at meter wavelengths. This 
synchrotron-type component has brightness temperature 



which is consistent with the integrated contribution from discrete radio sources (see Longair 1995 
for a recent review). The anisotropic galactic contribution is about three times larger. Such 
emission will cause intensity fluctuations from beam to beam, either from discrete radio sources 
or from variable galactic brightness. Hence the difficulty is to be able to discriminate between 
the IGM around a QSO with a brightness temperature of 0.01 K against a small fluctuation in 
the galactic and extragalactic background. Clearly, to remove this foreground contamination, 
some differential mapping technique must be used. The line radiation from high-z gas might, e.g., 
be distinguished against continuum emitting sources owing to its sharp and distinctive spectral 
feature, i.e., by making measurements at two frequencies sufficiently close that the continuum 
emission from discrete sources does not vary significantly (Bebbington 1986). 

In the following sections we will illustrate by a few examples the sensitivity of the detailed 
pattern of sky brightness at redshifted 21-cm wavelengths to the unknown thermal and ionization 
history of the universe. 




(47) 
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5.1. Quasars in a Background of UV Radiation at z 8 

In the theories for the origin of cosmic structure that are most popular today, massive 
objects grow by hierarchical clustering - low mass perturbations collapse early, then merge into 
progressively larger systems. In such a scenario, one naturally expects the high redshift universe 
to contain small scale substructures, systems with virial temperature well below 10 K. The 
resultant virialized systems will remain as neutral gas clouds unless they can cool due to molecular 
hydrogen, in which case they may form Pop III stars. These stars will then ionize the gas in their 
vicinity. It is less clear, however, whether or not a sufficient number of ionizing photons will be 
able to penetrate all of the surrounding neutral gas in the cloud to reach the external intergalactic 
gas and ionize it. Nonetheless, the stars will still produce a background of UV continuum photons 
near the Lya frequency which can escape. 

It is therefore of interest to consider a case in which QSOs turn on in a universe where the 
Lya scattering rate is P a ~ P t h everywhere, i.e., where 

J a = Jth - \^ k ^ « (lO—ergscm-s^Hz^sr^) (1±0) , (48) 
2A a vw c z \ 9 / 



(cf. eqs. [2C] and f22j). The Lya photons will propagate into uncollapsed regions of the IGM 



where the kinetic temperature in the absence of heating will be 

T K (z) « 0.026 K(l + zf (49) 

(Couchman 1985), well below Tcbr because of adiabatic cooling during cosmic expansion. The 
Lya flux will couple the spin temperature to the kinetic temperature, and T$ will be pulled below 
Pcbr- (Note that, in the absence of a radiative coupling agent like Lya , this could only happen 
in very high-baryon-density models, where collisions would remain important even at modest 
redshifts, Scott & Rees 1990.) 

As shown in § 4.1, however, the same Lya photons which mix the hyperfine levels will also 
heat the gas through recoil. Integrating the energy equation, assuming the Lya field is turned on, 
or more precisely reaches the level of equation (fi£|), at redshift z a , and using equations and 
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@, we find 

T K {z) « (0.026K)(l + z) 2 + (160K)(l+z) 2 [(l + z)- 5 / 2 - (1 + Zq )- 5 / 2 ]. (50) 

The characteristic timescale for heating the medium above the CBR temperature via Lya resonant 
scattering at the thermalization rate is 

At heat = -— ^ — A 10 «10 5 yr, (51) 

about 20% of the Hubble time at z ~ 8. The heating rate corresponding to equation (fSj) is 
shown in Figure 2. The result is a finite interval of time during which Lya photons couple the 
spin temperature to the kinetic temperature of the IGM before heating the IGM above the CBR 
temperature. If Lya sources turned on at redshifts z a < 10, this interval would present a window 
of opportunity in redshift space near z ~ 8 that would enable a large fraction of intergalactic gas 
to be observable at ~ 160 MHz in absorption against the CBR. 

In such a scenario, 21-cm emission on tens of megaparsecs scale would still be produced 
by the relatively rare, X-ray heated regions associated with QSO sources. If we denote by Ts, w 
the spin temperature of the "warm" (Ts jW 3> 2~cbr) gas, and by Ts )C the spin temperature 
of the surrounding, more pervasive "cold" (Ts iC < ^cbr) medium, as given by equation (|50|), 
the differential antenna temperature observed between these two phases of the IGM will be 
proportional to the effective efficiency 

A, = T CBR f ^£ - f^V (52) 

\ + S,c +S,w J 

Here, we distinguish between the neutral fractions in the cold and warm regions, to allow the 
warm zone to refer either to the inside or the outside of the ionized bubble. Even though the 
H I in the cold region is absorbing the CBR, by using the flux through the cold gas as a fiducial, 
the signature of the cold gas will appear as emission from the warm regions surrounding the 
quasars. The effective efficiency will remain much larger than unity, An ~ Tcbr/?S iC ^ 4, only for 
< 10 7 yr (cf. equation []50]]), to drop to Arj ~ 2 after ~5x 10 7 yr. Note that Ai] ~ 4 corresponds 
to temperature fluctuations of about 0.04 K. This flux is well below the level of the Bebbington's 
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(1986) radio survey, but would be reachable with an instrument like the GMRT in several hours 
of integration. 

For illustrative purposes, consider now the sphere of influence associated with a quasar 
(■^a,47 = 1) a s = 1-5) which turned on at z = 8.7, nearly synchronized with a Lya background of 
sources turned on at z = 9, and radiating isotropically. At z = 8, e.g., after 5 x 10 7 yr the QSO will 
have ionized a bubble of size rj ~ 7.5 Mpc, and heated the surrounding H I within r c ~ 15Mpc. In 
Figure 3a we show Ar] as a function of radius at various epochs, including the effect of Lya heating 
on the cold phase. The peak value at the I-front initially exceeds unity, but declines as T$ )C rises. 
The quasar has created a shell of 21-cm emission surrounding the H n region with a thickness of 
Ar ~ 7.5 Mpc. The hydrogen column density through the shell's width is ~ 2 x 10 21 cnr 2 , and its 
total H I mass is ~ 2 x 10 16 M . The angular size on the sky corresponding to the diameter of the 
emitting region is 

go(l + zfr c 

26 = — 2 degrees. (53) 

The shell will have a characteristic radial double peaked profile, as the gas in its interior will be 
completely ionized, as shown in Figure 3b. The velocity difference across each peak, due to the 
Hubble flow, is Av = Hq(1 + z) 3 ^ 2 Ar ~ 10, 000 km s -1 . Between the inner edges of the two peaks, 
At; is twice as large, since Ar ~ rj at this time. At 150 MHz, we would then observe two signals 
(the front and back of the shell) with flux density 



(54) 



(cf. eqs. [|45| and p2| ), extending over and separated by about 5 MHz, nearly the entire bandwidth 



of the receiver. For a fixed velocity width, 8S increases linearly with the mass of the shell. 

Let us analyze in some detail the detectability of a typical 21-cm emitting shell with the 
GMRT. The 30 antennas of the GMRT will be configured as a Y of 16 dishes with a central 
compact core of 14 dishes. Because the GMRT functions as an interferometer, it is insensitive to 
angular scales exceeding ~ X/d, where d is the separation of the dishes. If we wish to measure 
~ 2 degree-scale structures at 150 MHz (A ~ 2 m), the projected separation of antennas is required 
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to be about d < 50 m. Since the GMRT antennas have a diameter of about 45 m, the GMRT 
interferometric system will not be sensitive to 2 degree-scale structures as the antennas will start 
shadowing each other for projected spacings less than 45 m. 

An alternative strategy is to measure the total power by incoherent addition of the outputs 
of the 30 GMRT antennas, for which the effective collecting area will be 30 x 0.5 that of one 
antenna. For observations at 150 MHz the theoretical rms thermal noise of the GMRT output in 
incoherent mode is ~ lmJy for the following parameters: T(sys) ~ 600 K, A(eff) ~ 1000 m 2 for 1 
45-m dish, a bandwidth of 10, 000 km s -1 (5 MHz), and 5-hr integration for each of the 2 frequency 
bands separated by about 5 MHz (Swarup &; Malik 1996). A comparison with equation (^34|) shows 
that ~ 100 hours of integration will be required to detect an individual shell at 5<7 rms . This 
time is comparable to that proposed for searches for emission from collapsing H I clouds at high 
redshift (Swarup et al. 1991; Subramanian &: Padmanabhan 1993; Swarup 1996). However, it is 
not clear whether the above sensitivity can be realized in practice because of the man-made noise 
and variations of the instrumental gain and the galactic background temperature over 2 bands 
separated by 5 MHz. Note that, although the differential antenna temperature will be larger at 
earlier epochs because of the larger brightness constrast (i.e., less time to raise the temperature of 
the cold phase), the smaller radius of the QSO sphere of influence more than cancels the increase 
in emission efficiency. The signal-to- noise ratio improves slightly with increasing bubble size. 

To calculate the expected number of emitting shells in a radio survey, we may estimate their 
effective covering factor on the sky as 



where nnso is the (proper) number density of quasars, d& is is the angular-diameter distance, 
is the field of view of the radio antenna (about 3.1 degrees at 150 MHz), and Av measures the 
radial dimension covered by the survey. If the latter matches the shell thickness Ar, e.g., for a 
bandwidth of 5 MHz, the covering factor may be rewritten as 



C 



^n QS o4(0 + e) 2 — 



(55) 




(56) 
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where q2i is the volume filling factor of the 21-cm emitting bubbles. As @ 2 Ar/6 2 r c ~ 5, it follows 
from equation ( |56D that a shell will be detected after scanning about 20 fields with the GMRT if 
the filling factor of warm bubbles is ~ 1%. 



5.2. Isolated Quasars 

We shall now analyze the model problem of an individual, isolated QSO turning on at epoch 
z > 5 in an otherwise unperturbed IGM. In this case the only radiative coupling agent is the 
continuum UV radiation of the quasar itself, redshifted to Lya . From equations (^) and (^9) , 
it is useful to define a "thermalization distance" from a QSO, denoted by r t h, at which P a = P^. 
Numerically, this yields 

r th «(6Mpc)(i±^ L%. (57) 

For the best case of T a 3> Tqbr, the efficiency for 21-cm emission may be expressed, using 
equations fll7|), fll9|) , (|20|), and (|46|) , as rj = (1 + r 2 /r| ll )'~ 1 - At r = rth, ^ = 1/2- Comparison with 
equation (^) shows that the thermalization distance is comparable to the size of the H n region. 
Since r > rj is required for detection, rj > 1/2 may be achieved only for a brief redshift interval 
after the quasar turns on. This interval is so small for many quasars that their H II regions will 
still be in the luminal expansion stage (cf. equation 0]), so that they will not be preceded by 
a warming front and will thus not be visible in 21-cm emission. The efficiency will therefore 
generally lie in the range < rj < 1/2. 

Larger shells offer a slight advantage for detectability. This is because the flux of a shell 
emitting 21-cm radiation depends not only on its efficiency but on its size as well. The frequency 
integrated flux is proportional to 

(fiAz>) cff = 2vr J dv J dSOrj. (58) 

Consider a measurement made at a frequency corresponding to a depth in the emitting zone 
around the quasar just outside the I-front, and with a bandwidth Av c matched to the light radius 
r c . Denoting by # t h the angle on the sky subtended by r t h, integrating equation (p^) out to r c , 
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using the expression above for 77, gives (OAz/) e g ~ 27r0jf h Az/ c [l — (rth/^c) arctan(r c /r t h)], when 
r c 3> rj. While the signal increases initially like the cube of the size of the emitting region, it does 

3 /2 

so only linearly at later times. Since the noise increases like r c , the signal-to-noise ratio increases 
asymptotically only like r X J 2 . Allowing for a light radius that subtends an angle 6 C on the sky of 
1°, the flux is then 



over a bandwidth of Ai/ C « 11 MHz [(1 + z)/9] _1/2 (6> c /l ), comparable to the rms noise. It would 
thus take a few hundred hours of integration time on the GMRT to detect the shell at 5cr rms . A 
quicker strategy may be to search for an excess variance across the sky over the rms noise that 
these shells would produce in the observed flux per resolution (beamwidthx bandwidth) element. 
A next generation facility like the Square Kilometer Array Interferometer (Braun 1996) could 
easily map warm IGM bubbles, effectively opening much of the universe to a direct study of the 
reheating epoch. 



The spectral appearance at 21-cm of non- uniform gas at high- z had been previously 
investigated by Hogan & Rees (1979) and Scott & Rees (1990). These earlier papers showed 
that the 21-cm emission from high-z H I may display angular structure, as well as structure in 
redshift space; radio observations could therefore in principle yield "tomographic" mapping of 
protoclusters or other large-scale inhomogeneities at high redshifts, where the gas had not yet been 
reionized. The present work discusses how similar techniques may be used to reveal "patchiness" 
in the 21-cm emission due to non-uniformities in the spin temperature, rather than in the density. 
Such effects could tell us how (as well as when) the primordial gas was reheated; if the reheating 
were due to ultraluminous but sparsely distributed sources (e.g., the first QSOs), the resultant 
patches would have a larger scale than any gravitationally-induced inhomogeneities, and would 
therefore be less difficult to detect. The continuum radiation, redshifted to Lya , from an early 
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6. Summary 



-36- 



generation of stars would couple the spin temperature to the low temperature of the IGM for a 
transitory interval before heating the gas above the CBR temperature. The result is a "window of 
opportunity" during which it may be possible to detect the IGM in absorption against the CBR, 
and in so doing possibly reveal the first epoch of star formation. 

We have assessed three preheating mechanisms at early times: soft X-ray radiation from 
quasar sources, soft X-rays from virialized 10 10 — 10 11 Mq halos, and resonant scattering of 
background Lya continuum photons from an early generation of stars, and argued that radiation 
sources at high redshifts may prevent collapsing structures that cannot cool efficiently from being 
detected in 21-cm absorption against the microwave background. We have shown that scattering 
of Lya continuum radiation can mix the hyperfine levels of intergalactic H I within ~ 10 Mpc of 
an individual QSO, and that this will produce an appreciable 21-cm emissivity from the IGM in 
the quasar neighborhood. While the X-ray-heated bubbles will still survive as fossils, once the 
quasar has died, the Lya photons will rapidly adiabatically redshift out of resonance so that no 
emission will occur without an external supply of photons. 

If reionization occurred at redshift 5 < z < 10, these QSO spheres of influence might be 
detectable at meter wavelengths by the GMRT. Because of the very large scales involved, there is 
no difficulty in achieving a sufficiently narrow bandwidth and an adequate angular resolution with 
the GMRT; the main limitations are the low flux levels and the radio background subtraction. 

It is possible that both the reionization and heating of the IGM were due to an early 
population of stars on subgalactic scales (rather than quasars) at redshift z > 20, well beyond 
the reach of the GMRT. Even if they reionized the IGM at an epoch accessible to the GMRT, 
however, these sources would be less thinly spread, and the characteristic structures in the IGM 
generated by them would be smaller and less detectable, so that an alternative observing strategy 
than the one discussed here would be required. 

We acknowledge helpful conversations on the subject of this paper with P. Palmer, L. Spitzer, 
M. Voit, W. Watson, and M. White. We thank the referee, G. Swarup, for valuable comments 
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which helped to clarify the detectability of a neutral IGM with the GMRT. P.M. and A.M. were 
supported in part by NSF grant PHY94-07194 to the ITP, University of California, Santa Barbara. 
A.M. is also grateful to the W. Gaertner fund for its support. M.R. thanks the Royal Society for 
support. 



A. Lya Line Flux in the Neutral Medium 

We evaluate the Lya flux produced in the neutral H I outside an ionization front from 
recombinations within the H n region itself. The volume emissivity is given by 

e v = x e x p n^a c ^phv<j){v), (Al) 

where a^f P is the effective recombination coefficient to the 2 2 P level (Pengelly 1964). Integrating 
from just outside the ionization front into the H n region, and averaging over 4n sr, we get 

(with x e = 1.2, x p = 1, and including an additional factor of 1/2 since Lya photons arise only 
from within the H n region). In evaluating the integral, we relate the distance r — 77 from the 
ionization front to the frequency shift y from line center (in units of the thermal width) due to the 
Hubble expansion, by r — 77 = by/H, where b is the Doppler parameter. Because emission at line 
center is quickly redshifted out of the core, only the surface of the H 11 region contributes to the 
Lya flux. The resulting intensity can be expressed as 

J„_ 21 (z) « 0.03 *5f (5f^f (iii) 4 ". (A3) 

The recombination photons escaping from an expanding H II region produce then only a small flux 
of Lya photons. 

For completeness, we shall also compute the contribution from the local production 
of Lya photons in the warm (mostly) neutral material surrounding the ionized gas. As 
Lya photons are redshifted out of the Doppler core by the Hubble expansion, the local 
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background intensity becomes J a = (n^hc/AirH) [x e x p a^ p + ^ e ^Hi7eH+ £Ex / '(nnhu a )], where 
7eH ~ (2-2 x cm 3 s _1 ) exp(— II.84/T4) (Osterbrock 1989) is the collisional excitation rate of 
Lya photons by electron impacts, and the last term accounts for Lya collisional excitations from 
the nonthermal distribution of electrons generated by the soft X-rays. The factor £ ~ 0.4 (Shull & 
van Steenberg 1985) represents the net fraction of photon energy converted to Lya photons. For a 
temperature of T4 = 0.1 and allowing (optimistically) for x e ~ x p ~ 0.2, radiative recombinations 
dominate, and 

J a -2i(z) ~ 0.5 x e x p h 5Q I I I — ^— J . (A4) 

Recombinations and collisional excitations in a warm, largely neutral medium are therefore a 
negligible source of Lya photons. 



B. Heating by Resonant Scattering of Lya Photons 



We derive the rate of transfer of energy, equation ( p0| ) , from the radiation field to the atoms 
through recoil using the equation of radiative transfer. Neglecting expansion and source terms, 
the rate of change of the number density of photons in a given frequency band n v for an isotropic 
radiation field was shown by Basko (1981) to be given in the diffusion approximation by 



1 dn v I d 
-n m Av D — 



Av D -^- + 2fn u j 



(Bl) 



c dt 2 dv 

where nni is the number density of hydrogen atoms in the ground state, / = hv a /m-Rcb expresses 
the effect of recoil, and only terms of order v/c have been retained. This is only correct as long as 
kTx <C hv a . Note that, in a steady state, n v oc exp[— h{y — v a )/kTK] near the resonant frequency, 
as was first demonstrated by Field (1959b). 

The rate of change of energy in the radiation field due to recoil is then 

roc Q 

u a = hv a h a = n m cAv D h / dv v— {a u fn u ). (B2) 

Jo ov 

Performing the integration by parts, noting that a v is sharply peaked at the resonant Lya frequency, 
using n u = (Air / c) J u / hv , and setting E a = —u a /nm, we recover equation (30) of the text. This 
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result is valid in the presence of expansion and source terms as well, as may be shown from the 
generalization of equation (Bl) by Rybicki & Deh" Antonio (1994). We note that, to order v/c, the 
flow of energy is one-directional, from the radiation field to the atoms. Treating the reverse flow 
requires the inclusion of terms of order (v/c) 2 (including relativistic effects like aberration). We 
have also neglected stimulated emission. These effects become particularly important when the 
matter approaches thermal equilibrium with the radiation field. In this limit, the radiation field 
relaxes to a blackbody spectrum, and the net transfer of energy ceases. In the situations relevant 
to this paper, these effects are negligible. 
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Fig. 1. — Numerical calculation of the IGM surrounding a QSO (assumes qo = 0.5, H = 
50kms _1 Mpc -1 , Oigm = 0.05). The QSO was turned on gradually at z = 6 with a specific 
luminosity of L u = 10 31 ergs s _1 Hz -1 (y /Vl) -1,5 . The effects of secondary electrons are significant 
for soft X-rays, and have been included, along with radiative cooling mechanisms and Compton 
cooling off the CBR. (a) Shows the temperature of the IGM due to soft X-ray heating by the QSO. 
At z = 5.5, the IGM is heated to a temperature above that of the cosmic background radiation out 
to a radius of 23 Mpc. The light front is at 25 Mpc. (b) Shows the H n (solid) and He n (dashed) 
fractions. The ionized region surrounding the QSO extends to 7.5 Mpc at this time. 

Fig. 2. — Heating rate per hydrogen atom of the IGM due to soft X-rays from collapsed galaxy 
halos, in a flat qo = 0.5 and Hq = 50kms _1 Mpc -1 CDM-dominated cosmology (solid). Also shown 
is the required heat input over a Hubble time to match the CBR temperature, 1.5Tcbr(z)/H(z) 
(dotted). For z < 7, the rate of X-ray heating from collapsed halos is sufficient to raise the 
temperature of the IGM above that of the CBR. At higher redshifts, the IGM is too dense and the 
hot halos that can produce the required penetrating soft X-rays are too few for efficient heating. At 
these epochs, only part of the IGM will be heated above the CBR temperature. The heating rate 
E t h due to Lya resonant scattering (dashed) is shown for a background radiation field sufficiently 
strong to couple the spin temperature to the kinetic temperature of the IGM (see text). 

Fig. 3. — (a) The 21-cm efficiency is shown as a function of distance from a QSO with L a ^ = 1 
and as = 1.5, turning on at z = 8.7. A Lyct radiation field of intensity J a ,-2i = 0.8[(1 + z)/9] 
was turned on at z = 9 and partially couples the spin temperature of the neutral hydrogen to the 
kinetic temperature. The effect of heating of the IGM by Lya collisions is included. The efficiency 
is shown at time intervals of 2.5, 5, and lOxlO 7 yr after QSO turn-on, corresponding to z = 8.3, 
8.0, and 7.4, respectively. Because the IGM temperature is less than that of the cosmic background 
radiation at early times, a large efficiency is achieved in an X-ray warmed shell surrounding the 
quasar, (b) A map of the differential antenna temperature (K) as a function of angle (degrees) 
and velocity shift (10 3 kms -1 ) (differenced from 21-cm redshifted to 153 MHz), resulting from two 
quasars turning on at z = 8.7, as in (a). One quasar is at a comoving position corresponding to 
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z = 7.4 and the other at z = 8.0. The quasars are separated by 2.5° on the sky. Such a map could 
be constructed by making a sequence of angular maps in a given frequency band, here taken to be 
2000 km s -1 wide, and stepping in frequency until the full emitting region surrounding the quasars 
is recovered. In practice, a bandwidth of 5000 — 10000 km s" 1 will give a more feasible integration 
time for establishing a positive detection. The emission is measured relative to a reference beam 
well-separated from the quasars, either by frequency or by angle. 
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